Tomography of entangled massive particles 
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We report on tomographic means to study the stability of a qubit register based on a string 
of trapped ions. In our experiment, two ions are held in a linear Paul trap and are entangled 
deterministically by laser pulses that couple their electronic and motional states. We reconstruct the 
density matrix using single qubit rotations and subsequent measurements with near-unity detection 
efficiency. This way, we characterize the created Bell states, the states into which they subsequently 
decay, and we derive their entanglement, applying different entanglement measures. 
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Quantum state tomography^ allows the estimation 
of an unknown quantum state that is available in many 
identical copies. It has been experimentally demon- 
strated for a variety of physical systems, among them 
the quantum state of a light modeQ, the vibrational 
state of a molecule 0] and a single ionfl, and the wave 
packets of atoms of an atomic beam[5|. Multi-particle 
states have been investigated in nuclear magnetic reso- 
nance experiments!^ as well as in experiments involving 
entangled photon pairsf^, P|- However, no experiment to 
date has completely reconstructed the quantum state of 
entangled massive particles 0, Ull . In the con- 
text of quantum information theorv|l3|. multi-particle 
entangled states are considered to be an essential re- 
source for processing information encoded in quantum 
states (qubits). Many protocols require the determinis- 
tic creation of entangled states and the preservation of 
these states over times much longer than their produc- 
tion time. For experiments aiming at entangling qubits, 
the coupling of the qubits to the environment needs to 
be well understood so that its effect can be minimized. 

In this article, we describe the deterministic creation 
of all four two-ion Bell states, i.e. ^± = (|10) ± \Ql))^/2 
and $± = (111) ± |00))a/2, where |a;ia;2), G {0,1} 
represents the combined state of the two qubits. By to- 
mographically reconstructing the two-ion density matrix, 
we fully characterize these states and determine their en- 
tanglement by different means. We also measure the time 
evolution of these states and determine the decay of their 
entanglement. 

In our experiments, a qubit is encoded in a super- 
position of internal states of a calcium ion. We use 
the 151/2,771 = —1/2) ground state and the metastable 
I-D5/2, m ~ —1/2) state (lifetime r ~ 1 s) to represent the 
qubit states |1) and |0), respectively. Two ^^Ca"*" ions are 
loaded into a linear Paul trap having vibrational frequen- 
cies of (27r)1.2 MHz in the axial and (27r)5 MHz in the 
transverse directions. After 10 ms of Doppler and side- 
band cooling the ions' breathing mode of axial vi- 
bration at cjb — (27r)\/3x 1.2 MHz is cooled to the ground 
state |Gb) ( 99% occupation). Thereafter, the qubits are 



initialized in the |11) state. For quantum state engineer- 
ing, we employ a narrowband Ti: Sapphire laser which is 
tightly focussed onto either one of the two ions. By excit- 
ing the Si/2 to D5/2 quadrupole transition near 729 nm 
(7 — 0.16 Hz), we prepare a single ion in a superposition 
of the |0) and |1) states. If the laser excites the transition 
on resonance ("carrier transition"), the ion's vibrational 
state is not affected whereas if the laser frequency is set 
to the transition's upper motional sideband ("blue side- 
band"), the electronic states become entangled with the 
motional states |0b) and |lb) of the breathing mode. We 
use an acousto-optical modulator to switch between car- 
rier and sideband transition frequency and to control the 
phase of the light field il^ |15||, E| ■ An electro-optical 
beam deflector switches the laser beam from one ion, over 
a distance of 5.3 /^m, to the other ion within 15 /is. Di- 
recting the beam which has a width of 2.5 ^m (FWHM 
at the focus) onto one ion, the intensity on the neigh- 
boring ion is suppressed by a factor of 2.5 • 10~^. By a 
sequence of laser pulses of appropriate length, frequency, 
and phase, the two ions are prepared in a Bell state as 
described below. For detection of the internal quantum 
states, we excite the S1/2 to P1/2 dipole transition near 
397 nm and monitor the fluorescence for 15 ms with an 
intensified CCD camera separately for each ion. Fluo- 
rescence indicates that the ion was in the S1/2 state, no 
fluorescence reveals that it was in the D^/2 state. By 
repeating the experimental cycle 200 times we find the 
average populations of all product basis states |00), |01), 
|10) and 111). 

We create a Bell state by applying laser pulses to ion 
1 and 2 on the blue sideband and the carrier. Using the 
Pauli spin matrices crxjCyjCzj^l and the operators h and 

that annihilate and create a phonon in the breathing 
mode, we denote single qubit carrier rotations of qubit a 

by 
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and rotations on the blue sideband of the vibrational 
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Now, we produce the Beh state ^± = 
(1 10) ± |01))/\/2 by the pulse sequence C/^_,_ — 
iT!^(7r,±7r/2)i?2(7r,7r/2)i?+(7r/2,-7r/2) appHed to the 
|11) state. The pulse Rf {tt/2, —Tr/2) entangles the mo- 
tional and the internal degrees of freedom, the next two 
pulses i?^(7r, ±7r/2)_R2(7r, 7r/2) map the motional degree 
of freedom onto the internal state of ion 2. Appending 
another 7r-pulse, J7$^ = i?2(7r, 0)U^±, produces the state 
<i>± up to a global phase. The pulse sequence takes less 
than 200 fj,s. 

To account for experimental imperfections, we describe 
the quantum state by a density matrix p. For its ex- 
perimental determination we expand p into a superpo- 
sition p = XiOi of mutually orthogonal Hermitian 
operators Oi, which form a basis and obey the equation 
tr {OiOj) = Adij^^. Then, the coefficients Xi are related 
to the expectation values of Oi by Xi = ti{pOi)/4. For a 
two-qubit system, a convenient set of operators is given 
by the 16 operators cr-^'' cr^^\ = 0, 1,2,3), where 

CTj-"'' runs through the set of Pauli matrices l,ax,<Jy,crz, 
of qubit a. 

The reconstruction of the density matrix p is now ac- 
complished by measuring the expectation values {cr^'^'^ 

aj ) p. A fluorescence measurement projects the quan- 
tum state into one of the states \x1X2 



X, e {0,1}. By 
repeatedly preparing and measuring the quantum state, 
we obtain the average population in states \x\x-i) from 
which we calculate the expectation values of 
and ai^^ ® a^z^ . To measure operators involving Oy, we 
apply a transformation U that maps the eigenvectors of 
<jy onto the eigenvectors of ctz, i.e. UayU~^ = az, where 
U = i?(7r/2, tt). Similarly, the operator ax is transformed 
into az by choosing U — i?(7r/2, 37r/2). Therefore, all ex- 
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pectation values can be determined by measuring a 
a'z'^ or ® 'J'z^ . To obtain all 16 expectation val- 
ues, nine different settings have to be used (see table QJ. 
Since a finite number of experiments allow only for an 



estimation of the expectation values (a-^'* ® '^)''')p^ ^he 
reconstructed matrix pp> is not guaranteed to be posi- 
tive semi-definite 19]. Therefore, we employ a maximum 
likelihood estimation of the density matrix^^ jio] , fol- 
lowin g th e procedure as suggested and implemented in 
refs. |2(J, |21j. We parametrize the density matrix by 
the elements of the Cholesky matrix T related to p by 
p — T'^T. Since the reconstructed density matrix pfi 
can have negative eigenvalues, we cannot use p^ as a 
starting point for the optimization routine. Instead, we 
use pp = PpRP/tr{PpRP) where P projects onto the 
subspace spanned by the eigenvectors of pn having non- 
negative eigenvalues. 
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TABLE I: Additional pulses used for reconstructing a two- 
qubit quantum state. All expectation values needed for quan- 
tum state reconstruction can be determined by applying ei- 
ther one of the nine different transformations given in the ta- 
ble below and subsequently measuring the qubit state. Only 
the expectation values of the operators displayed in bold face 
are used to reconstruct the quantum state. 
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i?(7r/2,37r/2) 
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For the pulse sequence that is designed to produce the 
state ^+ = (|10) + |01))/-\/2, we obtain the density matrix 

shown in Fig.^. The fidelity F of the reconstructed 
state is F = {^+\pq,^\-^+) = 0.91. To produce the state 

~ (1 10) — |01))/-\/2, we change the phase of the side- 
band TT-pulse by tt and experimentally obtain the density 
matrix shown in Fig. QJd. For $± = (|ll)±|00))/\/2, 
we find the density matrices depicted in Fig. ^ and d. 

Having reconstructed a Bell state's density ma- 
trix, we can now check that the two qubits are in- 
deed entangled. It has been shown that a mixed 
state p of two qubits is entangled if and only if its 
partial transpose has a negative eigenvalue [l^, 

1231 . Not surprisingly, the partial transpose of 
the density matrix pqi^ (Fig- W^) h^-s eigenvalues 
{-0.42(2), 0.40(2), 0.49(2), 0.53(3)} close to the values of 
a maximally entangled state {—0.5,0.5,0.5,0.5}. Errors 
in the determination of the density matrix elements and 
of quantities derived from them occur mainly as a conse- 
quence of quantum projection noise. Systematic effects 
like pulse length errors or addressing errors (coherent ex- 
citation of an ion by stray light) play a minor role. We 
estimate the magnitude of quantum projection noise by 
a bootstrapping technique |24| where the reconstructed 
density matrix serves for calculating probability distri- 
butions used in a Monte-Carlo simulation of our experi- 
ment. 

Among the different measures put forward to quan- 
tify the entanglement of mixed bipartite states, the en- 
tanglement of formation, E, has the virtue of being an- 
alytically calculable from the density matrix of a two- 
qubit system. For a separable state, E = 0, whereas 
i? = 1 for a maximally entangled state. Using the for- 
mula given by Wootters.25|, we find E{p\i,^) = 0.79(4). 
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FIG. 1: (a): Real and imaginary part of the density matrix 
p*_^ that approximates = (|10) + |01))/\/2- The mea- 
sured fidelity is F^_^^ = {*+|pvi,^ |*+) = 0.91. (b): Real and 
imaginary part of the density matrix that approximates 
= (1 10) - |01))/^2. The measured fidelity is Fvi,_ = 0.90. 
(c),(d): Density matrix elements of (c) and (d). 
Here, = 0.91 and Fi>_ = 0.88. 



For the density matrices that nominally correspond to the 
Bell states and <i>±, we measure = 0.75(5), 

Eip^^) = 0.76(4) and E{p^_) = 0.72(5). 

All of these Bell states violate a Clauser-Horne- 
Shimony-Holt inequality. For the state , for example, 
we obtain \ {A) \ = 2.52(6) > 2, where we have introduced 

the operator A — ai^"* (^^ (^^x-z + c^i^^ (X" cf^^+z + ® (^^x^z ~ 
(J^P ® with ax±z = [cTx ± o-^)/\/2. 
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FIG. 2: Time evolution of the state *+ = (|10) + |01))%/2. 
A magnetic field gradient causes the relative phase between 
the 1 10) and the |01) parts to evolve linearly in time. Deco- 
herence leads to a loss of fidelity, (a) Phase Pm{t) for which 
the overlap Fm = max^({4'^|p*^ (r)|4'^)) between p*^(r) 
and states of type — (|10) + exp(i/3)|01))\/2 is maxi- 
mized, (b) Fidelity F = (^/s^, (*)|p(T)|*,3i, (r))) of the mea- 
sured density matrix p(r) that nominally corresponds to the 
state */3i,(r) = (|10) -I- exp(j/3L)|01))\/2 (filled circles). The 
state $+ = (111) + |00))%/2 (t riangles) is sensitive to fluctu- 
ations of the laser frequency and the magnetic field. Thus, 
its decay occurs on a much shorter time scale. Error bars 
account for quantum projection noise as well as systematic 
errors in the reconstruction process. They are derived from a 
Monte- Carlo simulation. 



Once a Bell state has been produced, we monitor its 
evolution in time by waiting for a time t before doing 
state tomography. We expect the Bell states ^± to be 
immune against collective dephasing due to fluctuations 
of the magnetic field or the laser frequency How- 
ever, they will only be time-invariant if the energy sep- 
aration huj between the qubit states |0) and |1) is the 
same for both qubits. A magnetic field gradient that 
gives rise to different Zeeman shifts on qubits 1 and 2 
leads to a linear time evolution of the relative phase be- 
tween the |01) and the |10) component of the 'i'± states. 
This is indeed the case in our experiments. We calcu- 
late the maximum overlap Fm = niax/3((5'^|p^^ (r)!^'/?)) 
between the density matrix p*_|_(t) and the states 5'/? = 
l/'\/2(|10)-|-exp(i/3)|01)) as a function of time. The phase 
f3m (t) for which the maximum overlap Fm is obtained is 
drawn in Fig. |2Ji. The phase changes linearly with time 
according to Pl = with LUp = (27r)170 Hz, revealing 
the presence of a field gradient of dB/dz = 0.6 G/cm 
in the direction of the ion string that transforms a ^4. 
state into a state within 3 ms and vice versajl^. The 
decay of the Bell state into a mixed state leads to a slow 
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FIG. 3: Decay of the Bell state *+ = (|10) + |01))\/2 into a 
statistical mixture as a function of time, (a) Entanglement of 
formation and (b) smallest eigenvalue of the partial transpose 
of the density matrix. The insets show the magnitude of the 
density matrix elements measured after a time of 2 ms and 
8 ms. 



decline of the fidelity F = (^/3t|p*+(''')|^'/3i) over time 
which is shown in Fig. ^p. We find that the states 
decay into a statistical mixture within 5 ms whereas the 
<i>± states already decay after about 200 /xs to a fidelity 
of F — 0.75. That decay can also be seen by drawing the 
entanglement of formation and the smallest eigenvalue 
of the partial transpose as a function of time as shown 
in Fig. 121 In some measurements, wc found even longer 
lifetimes of the state approaching 20 ms. 

With the methods demonstrated above we have deter- 
ministically created all four Bell states. State tomog- 
raphy yielded the complete information about the two- 
qubit quantum states. In future, this tomographic pro- 
cedure will provide an efficient tool for evaluating quan- 
tum gates and the effect of decoherence in multi-qubit 
systems. We gratefully acknowledge support by the Eu- 
ropean Commission (QUEST and QGATES networks), 
by the Austrian Science Fund (FWF), and by the Insti- 
tut fiir Quantcninformation GmbH. 



[2: 
[3 
[4 

[6: 
[7: 

[8 
[9 

[10: 
[11 



[12: 

[is: 



[14 

[is: 
[16: 

[17 

[is: 

[19 
[20 

[21 

[22' 
[23 

[24 
[25 
[26 



D. T. Smithey, M. Beck, M. G. Raymer, A. Faridani, 
Phys. Rev. Lett. 70, 1244 (1993). 

T. J. Dunn, I. A. Walmsley, S. Mukamel, Phys. Rev. Lett. 
74, 884 (1995). 

D. Leibfried, D. M. Meekhof, B. E. King, C. Monroe, 
W. M. Itano, D. J. Wineland, Phys. Rev. Lett. 77, 4281 

(1996) . 

Ch. Kurtsiefer, T. Pfau, J. Mlynek, Nature 386, 150 

(1997) . 

I. L. Chuang, N. Gershenfeld, M. G. Kubinec, D. Leung, 

Proc. R. Soc. London A 454, 447 (1998). 

A. G. White, D. F. V. James, P. H. Eberhard, P. G. 

Kwiat, Phys. Rev. Lett 83, 3103-3107 (1999). 

A. G. White, D. F. V. James, W. J. Munro, P. G. Kwiat, 

Phys. Rev. A 65, 012301 (2001). 

E. Hagley, X. Maitre, G. Nogues, C. Wunderlich, 
M. Brune, J. M. Raimond, S. Haroche, Phys. Rev. Lett. 
79, 1 (1997). 

Q. A. Turchette, C. S. Wood, B. E. King, C. J. Myatt, 
D. Leibfried, W. M. Itano, C. Monroe, D. J. Wineland, 
Phys. Rev. Lett. 81, 3631 (1998). 

D. Kielpinski, V. Meyer, M. A. Rowe, C. A. Sackett, 
W. M. Itano, C. Monroe, D. J. Wineland, Science 291, 
1013 (2001). 

F. Schmidt-Kaler et al.. Nature 422, 408 (2003) 

M. A. Nielsen, I. L. Chuang, Quantum Computation and 
Quantum Information (Cambridge Univ. Press, Cam- 
bridge, 2000). 

Ch. Roos et al., Phys. Rev. Lett. 83, 4713 (1999). 

F. Schmidt-Kaler et al., J. Phys. B: At. Mol. Opt. Phys. 

36, 623 (2003). 

H. Haffner et al., Phys. Rev. Lett. 90, 143602 (2003). 

Note that jO) and jl) are eigenvectors of az corresponding 

to the eigenvalues +1 and —1, respectively. 

U. Fano, Rev. Mod. Phys. 29, 74 (1957). 

Z. Hradil, Phys. Rev. A 55, R1561 (1997). 

K. Banaszek, G. M. D'Ariano, M. G. A. Paris, M. F. 

Sacchi, Phys. Rev. A 61, 010304 (1999). 

D. F. V. James, P. G. Kwiat, W. J. Munro, A. G. White, 

Phys. Rev. A 64, 052312 (2001). 

A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

M. Horodecki, P. Horodecki, R. Horodecki, Phys. Lett. 
A 223, 1 (1996). 

B. Efron, R. Tibshirani, Stat. Science 1, 54 (1986) 
W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
The value of the field gradient was confirmed in an inde- 
pendent measurement. 



[1] K. Vogel, H. Risken, Phys. Rev. A 40, 2847 (1989). 



